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o 

The electrodynamics with a Chern-Simons term p^AvF^" violates Lorentz and CPT symmetries 
with a non- vanishing p M . For a fixed vector p M , in this paper we point out that the energy-momentum 
£0. tensor of this theory coupled to the gravity minimally is symmetric but not divergence free, which 

consequently makes the gravitational field equation inconsistent. To preserve the consistency, we 
introduce a Chern-Simons term in the gravity sector with the coefficient determined by the Lorentz 
and CPT violating term in the electromagnetic field. Further we study the phenomenologies of the 
£ — ' model on the cosmic microwave background radiation and the relic gravitational waves. 
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Both the Standard Model of particle physics and Einstein's general relativity are locally Lorentz and CPT invariant. 
Probing their violations is an important way to search for new physics and in recent years has attracted a lot of 
interests. Usually, studies are focused on the phenomenologies in matter and gravity sectors separately. However 
there is an issue related to the consistency of the theory which needs to be examined. In a class of models of 
Lorentz symmetry breaking, fixed preferred frames are introduced [l|, The violation effects are formulated by 
introducing operators like p^K^ in the Lagrangian of the matter with p M being a fixed vector and the matter 
current. The existence of the fixed preferred frames violates general coordinate covariance also. If the Lorentz 
invariance in the matter sector is broken in this way, it is impossible to get an energy-momentum tensor of the matter 
t-H ■ which is both symmetric and covariantly conserved. When gravity is included, this makes the Einstein's equation 
lO ' G^ v = R^ v — \ j2Rg^ v = —SirGT* 11 ' inconsistent. Because G^ u is symmetric between the indices /i and v, and its 
f — . divergence is identically zero as the result of the contracted Bianchi identity. In this paper we provide a solution to 
this problem by modifying the gravity simultaneously in a non-covariant way, so that the new terms in the gravity 
ON . sector balance the Lorentz violating effects in the matter sector. 

As a concrete example, in this paper we consider a matter sector where the electrodynamics is modified by a Chern- 
J> , Simons term p^A v F^ v . This term breaks Lorentz and CPT symmetries if p^ is treated as an external field and the 
related phenomenology has been studied extensively in the literature, such as testing CPT of photons in astrophysics 
[H and cosmology 0,H, H, 0> B S EH EH EH • A salient feature of this modified electrodynamics is the rotations of the 
polarizations of propagating photons. The rotation angle Ax depends on the external field. Due to this feature, a part 
of E type polarization will be rotated to B type polarization for photons. This will generate TB and EB correlations 
in the power spectra of the cosmic microwave background radiation (CMB). Hence we can use the CMB experiments 
to test the Lorentz and CPT violation in this model. With homogeneous and isotropic rotation angle, in Ref. two 
of us with Feng and Li did the simulations on the measurement of A\ with the CMBPol and PLANCK experiments. 
We pointed out that in such experiments the EB spectrum will be the most sensitive probe of such Lorentz and CPT 
violation. In two of us with Feng, Xia, and Chen first found that a nonzero rotation angle Ax = —6.0 ± 4.0 deg is 
mildly favored by the CMB polarization data from the three- year Wilkinson Microwave Anisotropy Probe (WMAP3) 
observations 0, EE EEE3, El and the January 2003 Antarctic flight of BOOMERanG [H, HO, El . With the newly 
released five year data [3], the WMAP group gives Ax = — 1.7 ± 2.1 deg, which when combined with the BoomeRang 
data is [H Ax = -2.6 ± 1.9 deg. The result by the QUaD collaboration is Ax = 0.55 ± 0.82 ± 0.5 deg \u\ and most 
recently improved to Ax = 0.64 ± 0.5 ± 0.5 deg 

As mentioned above, for a fixed p^, the electromagnetic Chern-Simons term is not invariant under the coordinate 
transformation. The energy-momentum tensor is the same as that of the Maxwell theory, but it is not covariantly 
conserved because the equation of motion is modified. This will make the Einstein's equation inconsistent. To preserve 
the consistency, we modify the Einstein's gravity simultaneously by introducing a gravitational Chern-Simons term. 
As a result, the left hand side of the gravitational equation is modified by a four dimensional Cotten tensor which 
also has non-zero divergence and will match the divergence of the electromagnetic field on the right hand side. We 
further in this paper study the phenomenologies of this model on the CMB and relic gravitational waves. 

Our paper is organized as follows: in section II, we firstly review briefly the Maxwell-Chern-Simons theory and 
point out the inconsistence when including the gravity. We then introduce a gravitational Chern-Simons term and 
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demonstrate how the theory becomes consistent; in section III, we study the effects of our model on CMB polarizations; 
in section IV, we analyze the late-time evolution of relic gravitational waves and our result show the potential signals 
at high-frequency regime; section V is the summary. 



II. CHERN-SIMONS MODIFIED ELECTRODYNAMICS AND GRAVITY 

The Maxwell-Chern-Simons theory is the Maxwell electrodynamics modified by a Chern-Simons term: 



= -\f^ + ^ V" , (1) 

where is the electromagnetic field tensor and F^ v — l/2e^ vpa F pa is its dual, p p — V M #i characterizes the preferred 
frame and has dimension [E\. We use the signature (+, — , — , — ) for the metric. In the second line of Eq. (fTJ), 
an integration by part is used. To characterize a fixed preferred frame, p p is assumed to be a constant vector in 
spacetime, so 9\ = p p x p + C. Here C is a constant, but it only contributes a surface term to the Lagrangian and can 
be set to zero. We will ignore the sources for the electromagnetic field in this paper. When considering the minimal 
coupling to gravity, from the Lagrangian fl} , we get the energy-momentum tensor 

T<r = -4|^, (2) 

where g = —detg^ v . It is this energy- momentum tensor provides the source to the gravity and appears in the right- 
hand side of the gravitational field equation. Because the Chern-Simons term is a topological term, it does not depend 
on the metric and has no contribution to the energy-momentum tensor. So Tp is the same as that of the Maxwell 
theory, 

Tp" = ^F aP F ap g^ - F^ a F v a . (3) 
However, the equation of motion is modified as 

V/'"' = p p F^ . (4) 
After making use of the equation above, we find the energy-momentum tensor is not covariantly conserved, 

= -\p v F m F» a . (5) 

Substitute the tensor in Eq. ([3]) into the Einstein equation, it becomes 

= -8TrG(T£ v + T£") , (6) 

where Tfif is the energy-momentum tensor of other matter. In this paper we assume there are no other Lorentz 
violations except the Chern-Simons term in the photon sector, so Tfif is symmetric and divergenceless. From Eq. ([5]), 
we see the Einstein equation is not consistent. 

To solve this problem we introduce a Lorentz violating term in the gravity sector simultaneously. There are many 
Lorentz violating modifications of gravity discussed in the literature. In this paper, we only consider the Chern-Simons 
modification proposed by Jackiw and Pi [22}. So the total Lagrangian is 

c = t^g {r+ t rR) - \ F ^ F " V + °i F ^ v + £ ™ > ( ? ) 

where RR = \^ va ^ R (y pa pR p allu is the gravitational Chern-Pontryagin density and R a pa p is the Riemann tensor. 
Similarly the parameter 62 — q^x 11 , where q p is a constant vector with dimension IE'' 1 ]. The variation of the action 
S = J d^Xy/—gC with respect to the metric gives the modified Einstein equation [22J, |23| 



(8) 
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where 

= -^[\7 a 9 2 (e^\7 a R^ + e^V a R^) + CT V ' T 9 2 (e^ R^ af3 + R^ a0 )} (9) 
is the four dimensional Cotten tensor, and its divergence gives 

= \q v RR ■ (10) 
o 

So, the divergence of Eq. (J8]) gives the constraint 

q u RR = 16TrGp"F fia Fi* a . (11) 

The normal Einstein field equation is a second order partial differential equation. When modified by the gravitational 
Chern-Simons term, the field equation is promoted to third order due to the four dimensional Cotten tensor shown 
in equation ©. This means the constraint (fTTj) will not make the whole system overdetermined. But the solutions 
will be different from those obtained in general relativity. For example, in Ref. the authors pointed out in the 
case FnaF^ = that the Schwarzschild solution is still existent but the Kerr solution is not. 

III. THE EFFECTS ON CMB ANISOTROPIES 

In the previous section, we considered in Eq. ([7]) the Chern-Simons modified electromagnetic field with fixed 
preferred frame and introduced the gravitational Chern-Simons term to modify the gravity sector simultaneously to 
preserve the consistency of the theory. In this section we will study the phenomenology of this model on CMB. For 
keeping the rotation invariance, we assume only the temporal components of and are not vanished. During 
inflation the density F iLV F ilv is diluted and can be set to zero. So the gravitational Chern-Simons term is constrained 
to be RR = 0. But at the linear order, as shown in Refs. 0, [22j], the gravity is modified and the produced tensor 
perturbations have different intensities for different helicities. The gravitational wave has two independent polarized 
components denoted by + and x . Here it is more convenient to use the right- and left-handed circular polarized 
components: 

h R =±(h + -ih*), 

h L = -L( h + + ih x). (12) 

The primordial gravitational waves generated during inflation freeze out after exiting the horizon. The power spectra 
for different handedness are defined as follows: 

(h R *Qn)h R (k 2 )) - P/?* s (ki - k 2 ) , (/i i *(k 1 )/ l i (k 2 )) = P^ 3 (k! - k 2 ) , (ft«*(ki)/i £ (k a )) = . (13) 

As mentioned above, due to the gravitational Chern-Simons term, the produced P R and P^ are not equal. The 
discrepancy depends on qq and the Hubble constant H ln during inflation and can be denoted by the small parameter 
e = -(7r/2)q H in 0: 

Ph = \Ph{l - e) , P£ = \Ph{l + e) , P,f + Pt = Ph , Ph - Pk = -ePh ■ (14) 

As is well known, the primordial scalar perturbations can generate the temperature and -E-mode polarization 
perturbations in CMB. The tensor perturbations can generate £?-mode perturbations besides T and E. In the 
absence of the Chern-Simons term the generated _B-modes are not correlated with T and E. As mentioned above, 
with the gravitational Chern-Simons term, the non-vanished TB and EB correlations would be produced at the last 
scattering surface, which we will explain in detail in the following. 

The polarization of electromagnetic field is described by the Stokes parameters /, Q, U and V. For CMB physics, 
the Stokes V is usually neglected because the Thomson scattering cannot produce net circular polarizations. The 
intensity / is invariant under coordinate transformations, but Q and U are not. The combinations Q±iU behave like 
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spin-2 variables under the rotation. Given a map of temperature and polarization, we can expand the perturbations 
in terms of spin-weighted harmonic function as below 



I m 

(Q±iU)(n) = ^a±2,i m ± 2 yim(") • 

I'm 

The expressions for the expansion coefficients are 

a T ,i m = J dOF4(n)T(n) 
a± 2 ,im = I dn± 2 Y l * m (n)(Q±iU)(h) . 
Instead of a,2,im and a-i.im-, it is convenient to introduce their linear combinations 



(15) 



aE,lm — — (02,2m + a~2,2m)/2 
<lB,!m = i{a-2,lm ~ a-2,im)/2. 



The power spectra are defined as 



( a *X'.l'm' a Xdm) — Q X X Sl'lS m ' 



(16) 



(17) 



(18) 



with the assumption of statistical isotropy, where X and X 1 stand for T, E and B. 

In real space, using the spin raising and lowering operators 6 and d , it is useful to introduce two scalar quantities 
E{n) and B{n) defined as [25|] 



E(n) = -i [3 2 (Q + iU)+6 2 (Q-iU) 



E 



B(n) 



lm 

i 
2 



(Z + 2)! 



1/2 



Yi m (n) 



(1-2)1 

\3 2 {Q + iU)-d 2 (Q-iU) 



E 



(1 + 2)1 



(1-2)1 



1/2 



Y lm (n) 



(19) 



For each Fourier component, we can simply work in the coordinate frame in which fe || z and then integrate over 
all the Fourier modes. The generated temperature and polarization perturbations by gravitational waves can be 
expressed as [25j : 



A^ T) (T ,n,fc) = [(1 - ,i 2 )e 2 ^(fc) + (1 - ^ 2 )e- 2l<t, h L (k)] A^fa, 
(A( T) ± i A[ / T) )(r ,n,fc) = [(1 T n) 2 e 2 ^h R {k) + (1 ± tfe^h L (k)] a£\t , M , fe) , 



(20) 



where the superscript T denotes tensor perturbations and A are the obtained Polnarev variables [26j | by integrating 
the Boltzmann equations. We used the conformal time r and tq to denote present value, x = fe(ro — t) and fx = k ■ n. 



The quantities A^ ; , A k - and A v ^ ' are given by 



ST) 



A ( P(T ,n,k) = fi 2 )e 2 ^h R {k) + {\- fi 2 )e- 2 ^h L {k)]K { p(T ^,k) 



Af( To ,n,k) = 



(1 - fi 2 )e 2 ^h R (k) + (1 - ii 2 )e- 2l *h L {k)\ £(x)A ( p (r , /*, fc) 
(1 - ii 2 )e 2 ^h R {k) - (1 - M 2 )e- 2 ^ L (fc)l B^A^^rc/x, fc) 



(21) 



where the operators are defined as £(x) = —12 + x 2 [l — d 2 ] — 8xd x and B(x) — 8a; + 2x 2 9 x . From these equations one 
can show that the spectra TT, EE, BB and TE only depend on the sum of the primordial spectra of gravitational 
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waves, Ph = + Pfr. But the cross correlations TB and EB relies on the difference P? — P^ = —eP^, 

C XX{T) = ^ 2 I k 2 dkP h (k)\A { ^(k)] 2 



Cj E{T) = (4tt) 2 J k 2 dkP h (k)Ag\k)Ag\k) 
Cj B(T) = -(Anfe J k 2 dkP h (k)A ( ^ (fc)A^ (k) 

C EB[T) = -(4^) 2 e J k 2 dkP h (k)Ag\k)Ag\k) . (22) 

The transfer functions A 2? (k) are defined in Eq. (30) of Ref . [25[. We can see that the effect of gravitational Chcrn- 
Simons term on CMB is only to generate the TB and EB correlations and leave other power spectra unmodified. 
This result is consistent with that of Ref. [27| • 

The scalar perturbations also have contributions to TT, EE and TE. So, the total power spectra of CMB should 

be 

n TT _ n TT(S) . n TT(T) n EE _ n EE(S) , n EE(T) n TE _ n TE(S) , n TE{T) 

qBB _ qBB(T) c tb _ (jTB{T) C EB _ qEB(T) ^ 

After the perturbations of CMB generated at the last scattering surface, the polarizations of CMB photons are rotated 
under the influence of the electromagnetic Chern-Simons term. The rotation angle is Ax = (po/2)(to — Tjs) [H, where 
Tis is the conformal time of the last scattering. In this case, except TT, all other spectra are changed 



n TT,obs 
°i 


= cr 








n TE,obs 
°Z 


= Q T£ cos(2A x )- 


Cj B sin(2A X ) 






n TB,obs 
U l 


= Q T£ sin(2A x ) + 


Cj B cos(2A X ) 






n EE,obs 
°Z 


= Cf £ cos 2 (2A x )- 


fCf B sin 2 (2A X )- 


Cf B sin 


(4A X ) 


n BB,obs 
°i 


= Cf £ sin 2 (2A X )- 


- Cf B cos 2 (2A X ) + 


Cf B sin 


(4A X ) 


n EB,obs 


= isin(4A X )(Cf £ 


-C BB ) + C BB cos 


(4Ax) . 





(24) 

In the formulas above, the quantities with the superscript obs are those observed after the rotation. We can see that, 
in the standard case without the Lorentz and CPT violations, only TT, TE, EE and BB in the right hand side of 
Eqs. (|24p are produced. Consider the modifications, the gravitational Chern-Simons produce non-vanished TB and 
EB in the right hand side of Eqs. which are calculated by Eqs. (|22l) . Then the electromagnetic Chern-Simons 

term further rotated the E mode to B mode polarization and change all but TT spectra as showed in the left hand 
side of Eqs. ([24]). 



IV. ENERGY SPECTRUM OF GWB AT LATE-TIME EVOLUTION 



In this section we focus on the dynamics of gravitational waves at late-time evolution. As shown in above section, the 
Chern-Simons term does not contribute on the total primordial power spectrum of tensor fluctuations but only affect 
their propagations. Correspondingly, we expect the energy spectrum of GWB would be modified by the Chern-Simons 
term. 

We start by giving the equation of motion for the tensor fluctuations in Fourier space [24j. 

(1 - X s q -)h s k " + 2H(l - lx s q Q -)K' + (1 - X s q -)k 2 h s k = 16^Ga 2 ^ , (25) 

where the prime denotes the derivative with respect to the conformal time r = J T ^4|y. The subscript s represents 

the two polarizations, with X R = 1 and X L = —1. a s is the anisotropic part of the stress tensor, constructed by the 
spatial components of the perturbed energy-momentum tensor. We would like to neglect it first, and then consider 
its contribution on transfer function later. 
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A. Canonical representation 

To simplify the equation of motion for tensor fluctuations, we introduce a quantity v s k 28], which is given by 



^(r)=ayi-A^- (26) 
Then we redefine the variables of fluctuations as 

< = Hh% , (27) 
which are often viewed as generalized Mukhanov-Sasaki variables. The equation of motion for u| is given by 

< + (fc 2 -^VK=0- (28) 

K 

S I! 

Eq. (|28|) has an asymptotic solution when we neglect the last term -A- which implies \kr\ ^> 1, and it is strongly 

k 

oscillating like trigonometric functions. This feature coincides with an adiabatic condition, which corresponds to the 
case that the effective physical wavelength is deep inside the Hubble radius. Therefore, the modes can be regarded 
as adiabatic when they are staying in the sub-Hubble regime with |fcr| 3> 1, and we may impose a suitable initial 
condition in virtue of WKB approximation, 

' e~ lkT , (29) 



. 2k 

for cosmological fluctuations. 

Once we have resolved the solutions to the above equations, we can obtain the tensor power spectrum for the 
polarization mode h s k . The GWB we may observe today should be characterized by the energy spectrum, defined by 

f d(0\p GW (r)\0) 

yi GW (k,T) = —— — — , (30) 

Pc\t) a ink 

where Pgw{t) indicates the energy density of gravitational waves, and the parameter p c (r) is the critical density of 
the universe. Since the GWB we observed has already re-entered the horizon, its mode should oscillate in the form of 
sinusoidal function. Accordingly, we can deduce the relation between the power spectrum and the energy spectrum 
at the scales of interests as follows 

W*,r)= i i^ 5B X:W,r), (31) 
where we have used the Friedmann equation H 2 (t) = ^p-p c (r). 



B. Transfer function 



Now we analyze the evolution of tensor perturbations in the GWB nowadays. Since the primordial gravitational 
waves are distributed in every frequency, once the effective co-moving wave number is less than aH, the corresponding 
mode of gravitational waves would escape the horizon and be frozen until it re-enters the horizon. The relation between 
the time when tensor perturbations leave the horizon and the time when they return is a ou tHout = o-inHin- Therefore, 
we have the conclusion that, the earlier the perturbations escape the horizon, the later they re-enter it. Moreover, 
once the effective co-moving wave number is larger than aH, the perturbations begin to oscillate like the plane wave. 
In the following, we will establish the relation to relate the power spectrum observed today to the primordial one. It 
can be described by the transfer functions which are defined as follows, 

P h s (fc,r)=T s (fc,T)P^(fc,rO , (32) 

where Tj indicates the end of primordial inflation. 
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In order to make clear every possible ingredient affecting the cvolvcmcnt of the GWB, it is suitable and reasonable 
to decompose the transfer function into three parts as follows [29j. 



T s (k,T) = F°F°F° 

I H(t) |2 X(t) |2 X( T ) |2 



(33) 



In the above formula, hl(r) is the exact solution of Eq. (|25[k h k (r) is an approximate solution of Eq. (|25[) by 
neglecting the anisotropic stress tensor; and h k (r) is also an approximate solution which is equal to h s k {Ti) if fc < aH 
while equal to plane wave if fc > aH. 

First, from Eq. (|29[) one can see that after horizon re-entering, gravitational waves begin to oscillate with a decaying 
amplitude proportional to l/^(r). Therefore, from the definition of h s k we get 



^cos[fc(r-r fc ) + s ], k>aH 
h s k (ri), k <aH 



(34) 



where 4n depends on the initial condition, A k is the maximum of the amplitude of oscillations, and Tk is the conformal 
time when fc = aH. Since we require this function to be continuous, there must be a matching relation that /i|(Tj) = 
[A k cos(j) s ]/v k {Tk)- Based on these relations one can get the first factor Ff as follows 



Ft 



K( z k) 



cos 2 [fc(r - r k ) + (j) 8 ]/ cos 2 (j) s 



(35) 



where we introduce the redshift 1 + z = do/ a(r) in aim of showing the effects of the suppression as a result of redshift. 
The index "0" indicates today, and Zk is the redshift when the modes re-entered the horizon fc = aH. Note that the 
relation of Zk and fc can be given by the following equation 



^ 0) (l + z fc )exp 



1 + 2 . 



(36) 



where the sum over i includes all components in the universe. Since the contributions from dark energy and the 
fluctuations in radiation are very small, here we ignore them and then get 



1 + z k 



1 



*eq 




4(fc/fc ) 2 

(1 H~ Zeq)Qrn 



(37) 



ilm /^r- ' ■ The factor F{ describes the redshift-suppressing effect on the primordial gravitational 



where z eq = — 1 

waves. Since this factor shows strongly oscillating behaviour which is inconspicuous to be observed in the GWB, we 
usually average the term cos 2 [fc(r — 77-) + cj) s ] and obtain | for it. 

Second, when considering the influence of the background equation of state of universe on the re-entry of horizon, 
we focus on analyzing the factor F£. Since the background equation of state w varies very slowly, it is profitable to 



assume that the evolution of the scale factor is of form a — a^{^- 
tensor er s and the Chern-Simons modifications, we resolve Eq 



with a = 



3 again and then have 



To ignore the anisotropic stress 



J a _i(-fcr) , 



(38) 



where T is the Gamma function and J v is the i^-th Bessel function. If Ifcrl ^> 1, there is such a relation that 



ui (fc,r) |2 



T 2 («+^) / try 

7r y 2 1 



U\ (fc,Ti) I 

51, considering that the phase should be continuous, hence we 



— • cos 2 (fcr+3 L ). To match with Eq 
have the solution that when tensor fluctuations re-enter the horizon the conformal time Tk = — fr- A similar relation 
was obtained in [30I ] during primordial stage. Finally, the second factor F% is given by 



FS 



T 2 {a 



2n 

— I COS " 



(39) 
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The second factor shows that, when the gravitational waves re-enter the horizon, there is a "wall" lying on the horizon 
which affects the tensor power spectrum. 

Third, during the evolution of tensor perturbations, the nonzero anisotropic stress tensor a s would more or less bring 
some effects on the GWB. This effect is pointed out by Steven Weinberg[3lJ, and usually the primary ingredients 
are the freely streaming neutrinos which damp the amplitude of the tensor power spectrum. This damping effect 
just makes power spectrum of tensor fluctuation times a constant but do not change the dynamics of the GWB's 
evolution. A combination of analytic and numerical calculations performed in Refs. [32|, [H, [3J, [HI, [H[ suggests that 
F| = 0.803f 3 for the frequency of relic gravitational waves among 10 -16 i/z and 10 _10 Hz is in high precision. 

Eventually, we have discussed three kinds of leading corrections in the transfer function which make contributions in 
the evolution of the GWB. One can see that the modifications brought by the Chern-Simons term contribute mostly 
to the first factor F±. This is because the Chern-Simons term of gravity sector mainly affect the physics at high 
energy scale. 



C. Analysis of today's GWB 

Using the transfer functions, we are able to connect the primordial gravitational waves with what we observe today. 
To substitute Eqs. ([55]) . (|3"5|) , and the damping factor F3 into (f3"2")l . and making use of the primordial tensor spectrum 
(fT4|) we can give today's tensor power spectrum as follows, 



P^,T ) = T S {k,T )P S h {k,Ti 



1- X s q k(l + z k ) f* r 2 (a + i) (2^ 2a 



l-\ s q k (1 + Zk) 2 2tt V" 



x P* . (40) 



When the frequency of GWB is small enough, the time for the corresponding mode re-entering the horizon is close 
to today, and one can see that the tensor perturbation power spectrum would agree with the standard theory very 
well. Therefore, we expect to find signals in the high-frequency region. To make a comparison with the normal 
energy spectrum of relic gravitation waves, we assume the slow-roll parameter approaches zero during inflation and 
the potential of the inflaton is Vi n f ~ M 4 with the scale M = 5 x 10 15 GeV. Consequently we can obtain the 
semi-analytical form of the present energy spectrum of tensor perturbations with Chern-Simons modifications as 
follows, 

where the numerical calculations show / = 3.104 75 x 10 32 , and g equals to 2.15691 x 10 14 when the frequency is 
between 10" 16 and lO -10 ^ but takes the value 3.34395 x 10 14 outside this region. In the numerical computation, 
we have taken w = — 1 and ao = 1. We provide the numerical results in Fig. [TJ 

From Fig. \T\ one can see that the energy spectra of relic gravitational waves with Chern-Simons term coincide with 
that in standard gravity theory at low-frequency regime, but are suppressed at high-frequency regime. Moreover, 
the suppressions of energy spectra strongly depend on the value of the parameter qq. A similar suppression effect on 
energy spectrum of GWB due to a spacetime non-commutativity was found in [33] , but it takes place at low- frequency 
regime. We can understand this effect as follows. A Lorentz- violating term often brings an effective mass for gravitons 
which could suppress the energy spectrum. Therefore, this suppression depends on the energy scale of the effective 
mass term. Although for both non-commutativity and Chern-Simons modifications, the Lorentz violations take place 
at high-energy scales, the effective mass caused by non-commutativity appears at low-frequency regime but that 
brought by Chern-Simons term happens at high-frequency regime. 



V. SUMMARY AND CONCLUSIONS 



In this paper, we pointed out that the Lorentz violations induced by fixed preferred frames in the matter sector will 
make the gravitational field equation inconsistent because the energy-momentum tensor cannot be both symmetric 
and covariantly conserved. We provided a solution to this problem by modifying the gravity simultaneously. In a 
concrete model, we considered the Chern-Simons modified electrodynamics which breaks Lorentz and CPT invariance. 
Simultaneously the gravity is also modified by a Chern-Simons term. The phenomenologies of this model on CMB 
and the late-time dynamics of relic gravitational waves have been studied. For these modifications, the gravitational 
Chern-Simons term generates non-vanished TB and EB cross-correlations of CMB at the last scattering surface 
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FIG. 1: The plot of energy spectra of relic gravitational waves. The black solid line represents the energy spectrum without 
Chern-Simons modification; the red dash line gives the curve of energy spectrum with qo = W^'Mp 1 ; the green dot curve shows 
the energy spectrum with go = 10 5 Mp 1 ; and the blue dot-dash line shows the energy spectrum with qo = 10 Mp 1 . Here the 
potential of inflation is taken to be Vi n f — M 4 in which M = 5 x 10 15 GeV and M p is the reduced Planck mass. 

and leave others unchanged. After that, the electromagnetic Chern-Simons term will rotate the generated _E-mode 
polarizations to B-mode ones and change all but TT spectra. For the late-time evolution of the relic gravitational 
waves, we found that the Chern-Simons term mainly contribute to the amplitude of GWB inside the horizon. In 
this case, the energy spectrum of GWB is suppressed at high-frequency regime, which depends on an effective mass 
brought by the Chern-Simons term. From the current result, we notice that this effect is hard to be detected by 
recent experiments. However, the studies on gravitational models with various Lorentz violations [38l l39l | and their 
gravitational perturbations [13, EH are particularly of theoretical interests. This model deserves further studies. It 
will be interesting to seek for Kerr-Newmann type solutions for the spacetime metric in the region surrounding a 
charged and rotating celestial body. 

In the literature, there is another class of Lorentz violating theories in which the Lorentz symmetry is broken 
spontaneously, for example, the external vector is replaced by the derivative of a dynamical scalar field, p M = V^. 
In this case, we should include the kinetic term and the potential of <f> in the action. During the evolution of the 
universe, the scalar field develops a non-zero and the Lorentz symmetry is broken in this background. Such 
Lorentz violation has been used to generate the baryon number asymmetry observed in our universe [43 . I43L ItI . |45[ . 
The scalar field may be the dynamical dark energy L 43] or the curvature scalar [tH . 145^ . In this case, there is no 
problem of consistency and it is not necessary to modify the gravity theory [46l . |47| . For example, in the Maxwell- 
Chern-Simons theory considered in Eq. (|T]). if p r = V^r/* and the kinetic term and potential of <f> are included, the total 
energy-momentum tensor Tp + Tf is symmetric and divergence free even though Tp and T^ v are not covariantly 
conserved individually. The divergence of Tp U is canceled by that of T^ v . An interesting effect in this case is that 
the rotation angle of the photon is dynamical and spacetime dependent Q . 
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